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Abstract 

In this paper the structure properties of asymmetrical nuclear matter has been calculated em- 
ploying AVis potential for different values of proton to neutron ratio. These calculations have been 
also made for the case of symmetrical nuclear matter with C/V14, AV14 and ^Vig potentials. In 
our calculations, we use the lowest order constrained variational (LOCV) method to compute the 
correlation function of the system. 
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I. INTRODUCTION 



The interpretation of many astrophysical phenomena depends on a profound understand- 
ing of different parts of physics. Nuclear physics has an important role in determining the 
energy and evolution of stellar matter. Most of calculations for asymmetrical nuclear mat- 
ter has a close relationship with astrophysics. These studies are also potentially useful for 
understanding the effective nucleon-nucleon interactions in dense asymmetrical nuclear mat- 
ter, an important ingredient in nuclear structure physics, heavy ion collision physics as well 
as compact star physics. Nuclear matter is defined as a hypothetical system of nucleons 
interacting without coulomb forces, with a fixed ratio of protons and neutrons, and can be 
supposed as an idealization of matter inside a large nucleus. The aim of a nuclear mat- 
ter theory is to match the known experimental bulk properties, such as the binding energy, 
equilibrium density, symmetry energy, incompressi bility, etc. , starting from the fundamental 
two-body interactions (Pandharipande & Wiringa 1979| ). 

A good many-body theory for nuclear matter can be useful for studying the details of 
nucleon-nucleon interactions. The observed phase shifts from scattering experiments plus 
the properties of the only bound two-nucleon system, the deuteron, aren't enough to obtain 
a unique nucleon-nucleon potential. Nuclear matter studies can help us understand better 
exactly how the properties of the matter are affected by different elements of a potential, 
and what sorts of features are required to produce the observed saturation. Nuclear matter 



studies may also indicate w^ 
(Pandharipande & Wiringa 



lethe r a potential model for nuclear forces is workable or not 



1979|). 



The starting point for a microscopic theory of finite nuclei is to solve the infinite matter 
problem. A solution of the infinite matter problem would also be the first step in obtaining 
the equation of state for dense matter, which is necessary in the study of neutron stars. 
At the end, it is simply a very interesting many-body problem in its own right. Methods 
developed for it should be helpf ul in other dense quantum fluids such as liquid helium 
(Pandharipande & Wiringa |l979| ). 

The starting point for any nuclear matter calculation is a two-b ody potential that models 
the nucleon-nucleon interaction (Pandharipande Sz Wiringa 1979|). The first nuclear matter 



calculations were done by Euler 



1937l |). Very little was known about the interaction of 



nucleons at that time (Pandharipande & Wiringa 1979|). At the same time Yukawa potential 



2 



was formulated as: 



^ = 7^, (1) 



where 7 is a constant and fx is defined as jfjj 
the mass of vr meson) and r is the relative distance between two nuc 



^ (C is the speed of light and i s 



eons (Cohen |l971| ; 



Wong 2004j ). Several years later, Gammel, Christian and Thaler ( 1957l |) introduced a 
potential of the form: 

V = Vc{r) + VT{r)Su. (2) 
In Eq. ([2]), Vc{r) is the central potential, ^(r) is the tensor potential and 

5'l2 = 3(0-1 ■ f)((72 ■ f) - CTi • (72 

is the usual tensor operator. Then the potential was allowed to depend at most linearly on 
the relative momentum p, and a spin-orbit term was added to it. 



V = Voir) + VTir)Su + VUr)L . S. 



(3) 



Where L is the relative angular momentum and S is the tota l spin of the nucleon pair. This 
was the form originally proposed by Wigner and Eisenbud ( 1941| ). 

In 1962 the two most widely used poten tials were introduced. Both abandoned the Wigner 



form. The Hamada and Johnston 



1962| ) model had the form. 



where 



V = Vc{r) + VTir)S^2 + VLs{r)L . S + VLLir)L,2, 



L,2 = [5lj + i(ri-(r2)]L'' - (L.Sf 



(4) 



and the Yale potential was defined as (Lassila et al. 1962| ). 



V = Voir) + VTir)S,2 + ^L5(r)L . S + y,(r)[(L.S)2 + L.S - 



(5) 



In 1968 another potential was introduced by Reid ( 1968l |). This potential has a central 
term, Vc{r), for uncoupled states (singlet and triplet with L = J) and for coupled state s 
(triplet with L = J ± 1) has the form of Eq. In 1974, Bethe and Johnston ([197^) 

introduced a potential that had the general form of the Reid potential. BJ potential has a 
very hard core in {S, T) = (0, 0), (1, 1) channels. 
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Generally the above potentials are limited to a few operators and don't fit the data 
for all the scattering channels very well. In many-body calculations of nuclei and nuclear 
matter, it is suitable to represent the two nucleon interaction as an operator (Lagaris & 



Pandharipande 



198l|): 



(6) 



p 

where V'^i^rij) are functions of the interparticle distance Vij, and Ofj are suitably chosen oper- 
ators. The nucleon-nucleon (NN) interaction scattering data uniquel y show the occurrence 
of terms belonging to the eight operators (Lagaris & Pandharipande 198l| ): 



1, ai.aj,Ti.Tj, {ai.aj){Ti.Tj), Sij, Sijiji.Tj), (L.S)ij, (L.S)ij(rj.rj 



(7) 



in the Vij. Many nuclear matte r calculations have been done with potential models 
(Lagaris & Pandharipande 19811] )• Th is pot ential has two different models. One of them 
is Reid-V a (Pa ndharipande & Wiringa |l979| ) and the other is BJ-II Vg (Pandharipande & 



Wiringa 



1979| ) model. There is also a Vq model. The Vi=i^^ terms are neglected in the Vq 



model. The HJ Vq model is obtained by neglecting the L.S and quad ratic spin-orbit terms 
in Hamada and Johnston potential ( Pandh aripande & Wiringa 1979| ). while the GT-5200 
potential (Pandharipande & Wiringa 1979| ) is itself of a Ve form. 

Another A^A^ interaction model is Vi2- In this model, in addition to the 8 operators of 
Eq. ([7]), there is four momentum-dependent terms: 



p=9-12 



L'^{ai.aj),L^{Ti.Tj),L'^{ai.aj){ri.Tj] 



(8) 



The Vi2 potential like the Vq model has two different forms, which are Reid-Vi2 and BJ-II 
Vi2 (Lagaris & Pandharipande |l98lj ). 

In 1981 a phenome nologi cally two-nucleon interaction potential was introduced by Lagaris 
and Pandharipande ( 1981| ). This potential was obtained by fitting the nucleon-nucleon 
phase shifts up to 425 MeV in S, P, D and F waves, and the deuteron properties. It has 
two additional terms other than the operators in Eqs. (3) and (4) and is called as V14 or 
Urbana V14 {UV14) potential. 



Of 



p=13,14 



{L.S)\{L.Sf{n.r,). 



In f7Vi4 model, the two nucleon interaction is written as: 

p=l,14 



(9) 



(10) 



where V^{rij) is the well known one-pion-exchange interaction, Vf{rij) is an intermediate 
range interaction and Vg{rij) is a purely phenomenological short-range interaction. 

Ther e is a lso another form of V14 potential which was proposed by Wiringa and collabo- 



rators 



1984| ). It is called Argonne V14 (v4Vi4) potential. It has the general form of UV14 



potential. The difference between AV^ and UV14 models are in how the functions V^(rjj), 
Vfi^Tij) and Vg{rij) are defined. 

Traditionally, A^A^ potentials are formed by fitting np data for T = states and either np 
or pp data for T = 1 states. Unfortunately, potential models which have been f itted only to 
the np data often give not a good description of the pp data (Stocks & Swart |l993l |). even 
after applying the essential correlations for the coulomb interaction. By the same token, 
potentials fit to pp data in T = 1 states give simply a mediocre description of np data. Sub- 
stantially, this problem is due to charge-independence breaking in the strong interaction. In 
the pr esent work we use an updated version of the Argonne potential, AVis model (Wiringa 
et al. 1995| ). that fits both pp and np data, as well as low energy nn scattering parameters 
and deuteron properties. This potential is written in an operator format that depends on 
the values of S, T and Tz of the A^A^ pair. AVis potential includes a charge- independent 
(CI) part that has 14 operator components (as in AV14 model) and a charge-independent 
breaking (CIB) part that has three charge- dependent (CD) and one charge-asymmetric 
(CA) operators. The four additional operators that break charge-independence are given by 



=15-18 



(11) 



where 



is the tensor operator. In between the operators of Eq. ( ITTl) . the first three represent 
charge- dependence while the last one represents charge-asymmetry. 

In this paper, we use the lowest-order constrained variational method (LOCV) to cal- 
culate the correlation function of the nuclear matter. Primarily, the techni que of LOCV 
was us ed to s tudy the bulk properties of quantal fluids (Owen et al. |l977l |: Modarres & 



Irvine 



1979a| ). The method was later extende d to ca lc ulate the symmetry coeffic ient fo r 



;he sem i-empirical mass formula (Howes et al. ' l978a| , 1979l | ; Mo darres fc Irvine 1979a| , 



et al. 



1979b|), the properties of beta-stable matter (Modarres & Irvine 



1978b|), the surface energies of quantal fluids (Howes et al. [1978b|) and the binding 



1979al . |1979bj : Howes 



energies of finite nuclei (Bishop et al. 



1978 



; Modarres [198^). The LOCV method was 



further extended for finite temperature calculation and it wa s very succes s fully a pplied to 
neutron, nuclear and asymmetrical nuclear matter (Modarres 1993| . 1995| . 19971 ]) in order 
to calculate different thermodynamic properties of these systems. Recently, LOCV calcula- 
tions have been done for the symmetric nucle ar ma tter with phenomenological two-nucleon 
interaction operators (Bordbar & Modarres 1997(1) and the asymmetrical nuclear matter 
with AVis potential (Bordbar & Modarres jl998 |). The incompressibility of hot asymmet- 
rical nucl ear m atter have been also investigated within an LOCV approach (Modarres Sz 



Bordbar 



1998|). Very rece ntly, so me nucleonic systems such as the spin polarized neutro n 



matter (Bordbar & Bigdeli j2007a| ). symmetric nu clear rn atter (Bordbar & Bigdeli (2007b|), 
asymmetrical nuclear matter (Bordbar & Bigdeli 2008a| ). and neutron star matter (Bord- 



bar & Bigdeli 



2008a| ) at zero temperature have been studied using LOCV method with the 



d. Th e thermodynamic properties 



2008b|), symmet ric nu clear matter 



realistic strong interaction in the absence of magnetic fie _ 
of the spin pola rized neutron matter (Bordbar & Bigdeli 5 
(Bigdeli et al. 2009| ). and asymmetrical nuclear matter (Bigdeli et al. 2010| ) have been 
also studied at finite temperature in absence of the magnetic field. These calculations have 
been extended in the presenc e of rn agnetic field for the spin polarized neutron matter at zero 
temperature (Bordbar et al. 2011| ). The LOCV method is a fully self-consistent formalism 
and it does not bring any free parameter into the calculation. It considers the normalization 
constraint to keep the higher order terms as small as possible. The functional minimization 
procedure represents an enormous computational simplification over unconstrained methods 
(i.e., to parameterize the short-range behavior of corr elation functions) that attempts to go 
beyond the lowest order (Bordbar & Modarres 1998l |). 

In the present work, we intend to calculate the structure function of asymmetrical nuclear 
matter using the LOCV method employing UVu, AVu and AVig potentials. So the plan of 
this article is as follows: The LOCV method is described in Sec. [Tll Section [1111 is devoted 
to a summary of the pair distribution function and the structure function. Our results and 
discussion are presented in Sec. HVl Finally, summary and conclusions are presented in sec. 

El 
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II. LOCV FORMALISM FOR ASYMMETRICAL NUCLEAR MATTER 



We consider a trial many-body wave function of the form 

* = (12) 

where $ is a slater determinant of plane waves of A independent nucleons, F is an A-hody 
correlation operator which will be replaced by a Jastrow form. i,e., 

F^SHm, (13) 

i>j 

and 5 is a symmetrizing operator. The cluster expansion of the energy functional is written 
as 

-^l + ^2 + ^3 + ---. (14) 

The one-body term Ei for an asymmetrical nuclear matter that consists of Z protons and 
N neutrons is 

e,= y: (15) 

Labels 1 and 2 are used instead of proton and neutron, respectively, and kf = (Sn'^pi)^ is 
the Fermi momentum of particle i {p — pi + P2)- 
The two-body energy E2 is 

E2 = ^j: < tmi2)\tj - jt > (16) 

ij 

and 

^^^^^ = "1^^^^^^^' ^^^^^^^ + /(12)V(12)/(12). (17) 
The two-body correlation operator /(12) is defined as follows: 

a,p=l 

a = {J,L,S, T, Tz} and the operators 0^{ij) are written as 

or-'=i,(^+^^i2),(^-^5i2), (19) 

where S'12 is the tensor operator. We choose p = 1 for uncoupled channels and p = 2, 3 for 
coupled channels. The two-body nucleon-nucleon interaction 1^(12) has the following form: 

18 

y(12) = ^y^(ri2)Of2, (20) 



where the 18 operators that are defined as before, are denot ed by the labels c, u, r, err, t, 
tr, Is, 1st, 12, 12a, 12t, l2aT, ls2, Is2t, T, aT, tT, and tz (Wiringa 1984| ). By using correlation 
operators in the form of Eq. (ITS]) and the two-nucleon potential fr om Eg . (120]) . we find the 
following equation for the two-body energy (Bordbar & Modarres 1998|): 



—{-] 



E (2^+1): 



JLSTT^ 



-1) 



L+S+T 



X 



2m 



dr 



+ ir({'^^ ~ + ~ 3K.)(4T - 3) + {Vt - 3Kt) 
X [T(6T| - 4)] + 2K.T,}a«'(A;^r) + \Vi2 - 31^2. 



+ (T42r-3X^2ar)(4T-3)lc«'(A;^r) 



E 

j=2,3 



+ (-6z + 14) - (z - 1)V^,,](4T - 3) + [Vt + Kt(-6z + U)VtT\ 

X [T(6T,2 - 4)] + 2Vr.T,]af{kFr) + [V^2 + Vn. + (V^2r + V^2.r) 

X (4T - ?>)]cf{kFr) + [Vis2 + V5,2.(4T - 3)]d»'(fc^r))/«' 
2m r 

+ + 2Vi2 - 2Vn. - 'iVis2 + [(V^,. - 2l^2r - 2Vn.r - ms2r) 

X (4T - 3)]6^(^^r)/f + ) - ^f\\l{k,r) 



where the coefficients a^^\x), etc., are defined as 



(21) 



al^) fx) 



S'ix) 



[X 



[X 



x^Il,tA^), 

x^[/3/j_i,T,(a;) +7/j+i,T,(a;)], 
x2[7/j_i,T,(x) + ^/j+i,T,(a;)], 
a^^[/?23^j-i,T,(a;) - /323/j+i,r.(a;)], 
x^uiIl,tA^)^ 

x^[V2Ij^i,tAx) + '^2Ij+i,tSx)], 

x'^[^2Ij-i,tAx) + A2/j+i,T,(a;)], 
x^[^3/j_i,T,(a;) + A3/j+i,T,(a;)], 



(22) 
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with 



/3i 

A2 



J+1 

2J+ 1 



L(L + 1) U2 

J{J^ + 2J+1) 

2J + 1 
J3 + 2 J2 + 2 J + 1 

2J + 1 
J(J^ + J+1) 
2J + 1 



7 = 



J 

2J + 1 
1) 



2J(J + 1) 



23 



2J- 
^3 

A, 



1^3 



2J+ 1 
J3 + 2 J2 + 3 J + 2 



2J + 1 



6 



J( + J + 2) 

2J+ 1 
_ J(J^ + J + 4) 
^~ 2J+1 
J^ + 2J^ + 5J + 4 
2J + 1 



(23) 



and 



^J,T2(a;) = J dqPTz{q)Jj{xq). 
Ptz{q) is written as [tiz or = — | (neutron) and +| (proton)], 



for i 



Prziq) = -vrmin (/cf^^, ^^2 



F2 
Z2 



(24) 



(25) 



for g < i 



and 



PT,(g) = o 



for g > i fcf^i + k^z2 



The J/ fx) are the famihar Bessel functions. 



Now, we can minimize the two-body energy, Eq. ( |2T]) . with respect to the v ariations in 



the funct ions f' , but sub ject to the normahzati on con straint (Owen et al. [19771]; Modarres 
& Irvine [1979 j . (l979bj; Bordbar & Modarres (l998j) 

1 



,E<^jl4,(12)-f(12)Nj>a=0, 
where in the case of asymmetrical nuclear matter the function hxz (x) is defined as 



(26) 



hTAr) 



2 1 k[r 



1 



= ±1 
= 



(27) 
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In terms of channel correlation functions we can write Eq. ( 126|) as follows: 



1 



\L+S+T 



(28) 



X dr[hl^{kFr)- ff{r)\af{kFr) 







As we will see later, the above constraint introduces a Lagrange multiplier A through which 
all of the correlation functions are coupled. From the minimization of the two-body cluster 
energy we get a set of coupled and uncoupled Euler-Lagrange differential equations. The 
Euler-Lagrange equations for uncoupled states are 



+ 



m r 



K-3K + (K-3K,)(4T-3) 



(29) 



+ {Vt - 3v;t)[t(6t| - 4)] + 2v;.r, + A 



m 



Vl2 - SVl2a 



+ {Vl2r-SVl2ar){4:T-3) 



'■_a 

,(1)2 



= 0, 



while the coupled equations are written as 



+ 



(2)" 



K + K + 2Vt -Vls + (K + + 2Vtr 



Vi,r)i4:T - 3) + (Vt + + 2VtT)[Ti6T^ - 4)] + 2K.T, + A 



m 



^2 + ^2. + (^2r + ^2..)(4r-3) 



X (4r-3) 



aa 



{2f 



+ 



2^(2)2 ^ I ^2 



m 
2¥ 



V,s2 + V, 



Is2t 



- 2V,o - 2V, 



n2a 



- + (V^sr - 21^2r " 2V^2.r " 3T4,2r)(4T - 3) 

(3) = n 

a ^1 



(2) (3) 3a 
da da 



(30) 



,(3)" 



a, 



+ 



K + K - ^Vt - 2Vls + (K + Kr - ^Vtr 

2Vi,r){^T - 3) + (Vt + Kt - 4y,T)[T(6T2 - 4)] + 2K.T, + A 
m r 



1^2 + V2. + (V2. + V2.r)(4T-3) 



.(3)2 



Vls2 + Vi 



Is2t 



X (4T-3) 



+ 



(3)2 ' 2 (2) 



2 fi/a 



m 
2F 



V^, - 2V^2 - 2Vi2a 



(31) 
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- 3V^,2 + {Vlsr - 2Vl2r ' 2Vl2ar " 3V;,2r)(4r - 3) 



X 



hi 



,(2) 



(2) (3)^" 



where 



9l-\kpr) = f!^\r)a'^\kFr). (32) 

The primes in the above equation means differentiation with respect to r. As we pointed out 
before, the Lagrange multipher A is associated with the normahzation constraint, Eq. (!28ll . 
The constraint is incorporated by solving the Euler-Lagrange equations only out to certain 
distances, until the logarithmic derivative of the correlation functions matches those of (?") 
and then we set the correlation functions equa l to hx^ir) (beyond these state-dependence 
healing distances) (Bordbar & Modarres 1998l |). Finally, by solving the above differential 
equations (Eqs. (129|) . (l30l) and (l3T]) ) numerically, we obtain the correlation functions. 



III. STRUCTURE FUNCTION 



There are two types of structure functions, dynamic S{^^,w), and static 5'(k) st ructu re 
functions. They measure the response of the system to density fluctuations (Feenberg 1969| ). 

T he sta tic structure function of a system consisting of A particles is defined as (Feen- 
berg (l969|): 

S(k) = 1 + ^/ rfVidV2e^^-'--pi(ri)pi(r2)[(/(ri,r2) - 1], (33) 

where pi(r) is the one-particle density and (/(ri,r2) is the pair distribution function. In 
infinite systems, pi(r) is constant (= p) and is a function of the interparticle distance 
^"12 = ki — r2|, therefore Eq. fl33l) takes the following form, 

5(k) = l + p/e"^-'--[(7(ri2)-l]rfVi2. 



For calculating the pair distribution 
expansion of (7(^12) as follows (Clark 



(34) 



unction, we use the lowest order term in the cluster 



1979|), 

g{.ri2) = f {r 12) gF{,r 12), 



(35) 



where /(ri2) is the two-body correlation function and (?f('"i2) is the two-body radial distri- 
bution function of the noninteracting Fermi-gas, 



gF{ri2) = 1 I'^ikpru)- 



V 



(36) 
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In the above equation, v is the degeneracy factor, and l{x) = 3x ^{sinx — xcosx) is the 
statistical correlation function or the slater factor. 

IV. RESULTS AND DISCUSSION 
A. Correlation function 

In Fig. [H we have plotted our result for the correlation function of symmetrical nuclear 
matter versus internucleon distance (ri2 = r) employing UVu, AV^ and AVis potentials 
at density p = 0.16 fm~^. Here the correlation functions are calculated from average over 
all states. We can see that the correlation function is zero at the internucleon distance 
r < 0.06 fm for the three potentials. This distance represents the famous hard core of the 
nucleon-nucleon potential. When the internucleon distance increases, the correlation also 
increases until approaches to unity, approximately at r > 3.8 fm. This means that at r 
greater than the above value, the nucleons are out of the range of nuclear force (correlation 
length). The value of correlation for AVis potential has a maximum greater than unity 
and then approaches to unity. However, for UVu and AVu potentials, there is no such 
a maximum. In Fig. [21 we have plotted the correlation function of asymmetrical nuclear 
matter employing AVis potential for different values of proton to neutron ratio {pnrat = 
0.2, 0.6, 1.0) at different isospin channels {nn, np, pp). From this Figure, it can be seen that 
for all values of pnrat, the correlation functions of nn and pp channels have the maximums 
greater than unity, whereas at np channel, there is no such a maximum. This means that at 
pp and nn channels, the nucleon-nucleon potential is more attractive than at np channel. We 
can see that at nn and pp channels, the maximum values of correlation function decrease by 
increasing pnrat. We have found that at pp and np channels, the correlation length decreases 
as pnrat increases, while at nn channel, by increasing pnrat, the correlation length increases. 
In addition, for each pnrat, the value of the correlation length at pp channel is greater than 
that of np channel, and the correlation length at nn channel has a greater value than pp 
channel. These have been clarified in Table [T] in which the values of the correlation length 
for different values of pnrat at different isospin channels have been presented. 



12 



B. Pair distribution function 



We know that the pair distribution function, g{r), represents the probabihty of finding 
two particles at the relative distance of r. In Fig. [31 we have plotted our results for the pair 
distribution function of symmetrical nuclear matter versus internucleon distance with UV14, 
AVu and AVis potentials at density p = 0.16 fm~^ . Our results are in a goo d agreement with 
those of others calculations employing the Reid potential (Modarres 1987|). Figure [3] shows 
that for r in the range of 1.1 fm to 3.4 /m, the pair distribution function corresponding 
to AVis potential is greater than those of UVu and AVu potentials. This is due to the 
behavior of two-body correlation as mentioned in the above discussions. In the Fermi gas 
model due to the absence of interaction between nucleons, the pair distribution function 
is not zero even in the small internucleon distances as shown in Fig. [31 But in the real 
system, in which there is interaction between nucleons, the value of g{r) at r < 0.06 fm is 
zero for the three potentials. The same as for the case of correlation function, this distance 
represents the hard core of the nuclear potential. From Fig. [3l it can be seen that the 
value of g{r) increases as the internucleon distance increases and finally approaches to unity, 
approximately at r > 4 fm. In Fig. [H we have plotted the pair distribution function 
of asymmetrical nuclear matter employing AVis potential at different values of proton to 
neutron ratio (pnrat) for p = 0.16 fm~^ and different isospin channels {nn, np, pp). We 
can see that at all channels, by increasing pnrat, the pair distribution function decreases, 
corresponding to decreasing of the correlation. Besides, from Fig. [H it can be seen that for 
each pnrat, the pair distribution functions of nn and pp channels have identical behaviors, 
while at np channel, g{r), behaves differently compared to the other two channels. These 
are corresponding to the behavior of correlation function at these channels. 



C. Structure function 



In Fig. [5l we have plotted our results for the structure function of symmetrical nuclear 
matter versus relative momentum (fc) with UV14, AV14 and AVis potentials at density p = 
0.16 fm~^. There is an overall agre ement between our results and those of others calculated 
with the Reid potential (Modarres 1987l |). From Fig. [5l it is seen that the nucleon-nucleon 
interaction leads to the reduction of the structure function of nuclear matter with respect 
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to that of the non-interacting Fermi gas system. In Fig. [6l we have plotted the structure 
function of asymmetrical nuclear matter with AVis potential at different isospin channels 
{nn, np, pp) for different values of proton to neutron ratio (pnrat) and p = 0.16 fm^^ . It 
is seen that similar to the pair distribution function, the structure function of nn channel 
is like that of the pp channel, especially at higher values of k. We have found that this 
similarity becomes more clear as pnrat increases. However, there is a substantial difference 
between structure functions of np channel and pp and nn channels. 

V. SUMMARY AND CONCLUSIONS 

Using the lowest order constrained variational (LOCV) method, we have computed the 
correlation function, the pair distribution function and the structure function of the symmet- 
rical and asymmetrical nuclear matter. In order to investigate the effect of nucleon-nucleon 
interaction on the properties of nuclear matter, we have also computed the pair distribution 
function and the structure function of noninteracting Fermi gas. Here, we have used AVig, 
potential to represent the nucleon-nucleon interaction for the asymmetrical nuclear matter. 
These calculations have been done at different isospin channels. In the case of symmetri- 
cal nuclear matter, the calculations have been done with UVi4^, AV14 and AVis potentials. 
There is an overall agreement between our results and those of others calculated with the 
Reid potential. It was seen that the nucleon-nucleon interaction leads to the reduction of 
the structure function of nuclear matter with respect to that of the non-interacting Fermi 
gas system. We have found that at np and pp channels, the correlation length decreases as 
the proton to neutron ratio (pnrat) increases, while at nn channel, by increasing pnrat, the 
correlation length increases. However, the behavior of the pair distribution function at np 
channel is considerably different pair from those of other two channels. This is due to the 
difference between the behavior of correlation functions of these channels. It was indicated 
that for higher k and pnrat, the structure functions of nn and pp channels are identical, 
corresponding to the similarity between the pair distribution functions of these channels. 
We have also shown that the structure function at np channel was different from those of 
nn and pp channels. 
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TABLE I: The correlation length of asymmetrical nuclear matter employing AVis potential for 
different values of proton to neutron ratio at different isospin channels {nn, pp and np). 



pnrat 


correlation length (/m) 


nn 


np 


pp 


0.2 


2.95 


2.09 


2.18 


0.6 


3.36 


1.97 


2.11 


1.0 


3.39 


1.94 


2.06 
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r(fin) 



FIG. 1: The correlation function of symmetrical nuclear matter employing UVu, AVu and AVis 
potentials. The correlation functions have been calculated from average over all states. 
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FIG. 2: The correlation function of asymmetrical nuclear matter employing AVi^ potential for 
p = 0.16 fm~^ and different values of pnrat at different isospin channels (nn, pp and np). 
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r(fm) 



FIG. 3: The pair distribution function for symmetrical nuclear matter calculated with UVi/^, AVu 
and AVi8 potentials at density p = 0.16 fm~^. The pair distribution function corresponding to 
the fermi gas is also brought for comparison. 
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FIG. 4: As Fig. [21 but for the pair distribution function of asymmetrical nuclear matter. 
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FIG. 5: The structure function of symmetrical nuclear matter with UV14, AVu and ^Vis potentials 
at density p = 0.16 fm~^. 
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FIG. 6: As Fig. [2l but for the structure function of asymmetrical nuclear matter. 
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